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ABSTRACT 


A  comparison  of  the  Surface  Integral  and  Adjoint  Difference  techniques 
are  presented  in  this  report.  This  comparison  is  made  from  a  brief  theore¬ 
tical  development  of  each  method  and  by  the  results  obtained  using  each 
method  for  several  radiation  transport  problems  of  the  clas3  of  a  target 
located  in  air  (or  air-ground  interface)  far  removed  from  a  point  neutron 
source.  Mont  targets  considered  thus  far  do  not  contain  fissile  materials; 
however,  it  is  demonstrated  that  the  calculational  tools  employed  in 
this  study  are  applicable  to  such  targets. 

The  theoretical  considerations  lead  to  the  conclusion  that  an  error 
is  associated  with  the  Surface  Integral  technique,  but  this  error  will 
probably  be  small  for  the  class  of  problems  considered  herein .  In  contrast 
to  this,  the  Adjoint  Difference  technique  is  "exact".  The  numerical 
results  presented  agree  with  these  predictions. 
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X.  INTRODUCTION 


The  objective  of  this  project  was  to  develop  efficient  calculations! 
techniques  to  predict  the  absorbed  radiation  within  geometrically  complex 
targets,  e.g.  re-entry  vehicles.  The  targets  are  normally  far  removed 
from  the  source  responsible  for  the  radiation  field;  therefore  the  overall 
proolem  can  be  described  as  a  geometrically  complex  deep  penetration 
problem. 

In  this  study,  it  was  assumed  (for  convenience)  that  the  target  was 
located  in  air  or  at  an  air-ground  interface.  The  overall  problem 
was  broken  up  into  three  separate  steps: 

1.  The  calculation  of  the  free  field  radiation  (deep  penetration 
problem)  Jn  the  air  (or  air-ground  interface)  with  the  target 
absent. 

2.  The  calculations  associated  with  the  geometrically  complex 
target  without  explicitly  accounting  for  the  source. 

3.  The  combination  of  the  results  from  the  first  two  steps  to 
obtain  the  quantity  of  interest  for  a  particular  source- 
target  geometric  configuration. 

% 

The  first  step  can  be  described  as  a  deep  penetration  problem,  which  is 

geometrically  simple  (one  or  two-dimensional).  Accordingly,  this 

problem  was  solved' using  deterministic  techniques,  i.e.,  the  discrete 
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ordinate  codes  ANISN  and  DOT  .  The  second  step  is  a  geometrically 
complex  problem,  but  usually  not  one  of  deep  penetration.  Hencs,  rather 
straightforward  Monte  Carlo  techniques  are  well  suited  for  the 
calculations!  tool.  The  code  MQRSE^  was  used  here.  A  code  was  vritten 


to  implement  the  third  step  for  demonstrative  purposes.  This  implementation 
of  this  last  step  of  the  overall  procedure  becomes  somewhat  involved 
for  arbitrary  orientations  of  the  target  relative  to  the  source.  The 
complications  are  introduced  by  the  fact  that  the  reference  coordinate 
system  used  for  the  calculations  associated  with  step  one  are  different 
than  the  reference  coordinate  system  used  for  step  2. 

The  calculation  procedures  employed  require  the  solutions  to  the 
forward  neutron  transport  equation  for  step  one  and  the  adjoint  transport 
equation  for  step  2.  The  adjoint  transport  equition  was  solved  with 
two  different  sets  of  boundary  conditions  resulting  in  tvro  different 
ways  of  obtaining  the  quantity  of  interest,  i.e.,  the  Surface  Integral 
Approximation  technique  and  the  Adjoint  Difference  technique.  Each  of 
these  methods  will  be  presented  below. 

The  theoretical  development  and  the  associated  computational  procedures 
for  both  the  Surface  Integral  Approximation  and  Adjoint  Difference 

i 

techniques  are  reported  in  detail  in  a  doctoral  dissertation  which  was 
completed  in  connection  with  this  project  and  reported  in  Ref.  4. 
Furwhermore,  sufficient  examples  are  included  in  Ref.  1*  to  demonstrate 
the  applicability  of  each  technique.  This  report  will  include  only  the 
more  important  aspects  of  the  theoretical  results.  The  reader  interested 
in  the  more  detailed  development  and  implementation  is  referred  to  Ref 
4  (R-4),  For  completeness,  the  results  of  all  test  problems  which  have 
been  considered  will  be  included  in  this  report  (some  of  these  are  also 
presented  in  R-4) . 


II.  THEORETICAL  CONSIDERATIONS 

A .  Introduction 

In  this  section,  the  problem  of  computing  the  effect  of  interest 
within  a  target  removed  from  the  source  responsible  for  the  radiation 
field  is  considered.  The  formulation  is  presented  so  that  the  effect 
of  interest  can  be  obtained  by  the  three  distinct  steps  outlined  in 
the  preceeding  section  using  both  the  Surface  Integral  and  the  Adjoint 
Difference  techniques.  The  development  here  is  kept  brief  since  the 
more  detailed  development  is  presented  in  R-4. 

B.  Problem  Statement 

Let  A  be  the  effect  of  interest  and  R(p)  be  the  response  function 
(i.e.,  the  contribution  to  A  due  to  unit  angular  flux),  then  the  problem 
is  to  solve  the  integral  equation 

A  =  /^(p)R(p)dP  (1) 

where  ^(p)  is  the  angular  flux.  The  angular  flux  X3  the  solution  to 
the  Boltzmann  transport  equation 

Hl(p)<i>i(p)  =  8(5),  (2) 

* 

vV tore  S(fl)  is  the  source,  H^(p)  is  the  integro-differential  operator  of 
the  Boltzmann  equation  (see  R-4) ,  and  subscript  "x"  refers  to  the  medium 
with  the  target  located  in  air. 

C.  The  Surface  Integral  Approximation 

For  the  Surface  Integral  method,  a  closed  surface,  denoted  by  r  , 

O 

is  selected  aprlori  which  encloses  the  target  and  excludes  the  source. 


Then  the  effect  of  interest  can  he  obtained  by  solving  the  integral 
equation 

A  =  -  /  (S-n) Mr-r  )  $  (p)**(p)dp.  (3) 

Pi  1 

In  Eq.  3,  p..  represents  the  phase  space  interior  to  the  surface  r  ,  n 

-L  8 

is  the  outward  directed  normal  to  the  surface  rg,  the  subscript  "in" 
represents  inward  directed  0  only,  and  $**(p)  is  the  solution  to 

HJ  (p)  **(p)  «  R(p).  (1>) 

In  Eq.  It,  K#(p)  is  the  operator  adjoint  to  the  operator  H^p)  of  Eq.  2, 
and  the  boundary  condition  is 

<f>*(p)  -  0  (5) 

for  (n*f2)  >0  at  the  surface  r  , 

s 

There  are  no  approximations  associated  with  Eq.  3,  but  this  formu¬ 
lation  for  the  effect  of  interest  is  of  little  practical  value  since 
it  is  still  necessary  to  solve  Eq.  2  for  <^(p).  However,  we  need 
<j>^(p)  at  the  surface  rg  for  inward  directions  only. 

To  introduce  the  Surface  Integral  Approximation,  consider  the 
solution  to 

H2(p)  «J»2(p)  =  S(p),  (6) 

where  Hg(p)  is  the  integro-differential  operator  of  the  Boltzmann 
equation,  S(p)  is  the  same  source  as  in  Eq.  2,  and  subscript  2  refers 
to  the  medium  with  the  'rget  removed.  Assume  that  the  ^(p)  of  Eq.  3 
can  be  adequately  approximated  with  ^(p)*,  then 

*We  are  assuming  the  inward  directed  component  of  the  ilux  at  the 
surface  ?  is  uneffected  by  the  presence  of  the  vehicle. 


(7) 


.*  =  •-/  5(r-rg)  <j>2(p)  4>*(p)dp 

h 

which  is  the  Surface  Integral  Approximation. 

To  summarize,  the  Surface  Integral  Approximation  consists  of  the 
following  steps: 

1.  Calculation  of  the  free  field  radiation  in  the  medium  (air) 
with  the  target  absent,  i.e.,  solution  to  Eq.  6. 

2.  Calculation  of  ^(p)  (solution  to  Eq.  4)  vith  boundary 
conditions  given  by  Eq,  5* 

3.  Carry  out  the  indicated  surface  integration  of  Eq.  7* 

The  relative  location  of  the  source  to  the  target  is  accounted  for  in 
Step  35  therefore,  many  relative  locations  of  the  source  and  target 
can  be  analysized  by  a  single  completion  of  steps  1  and  2  and  repeated 
application  of  step  3. 


D.  The  Adjoint  Difference  Method 
For  the  adjoint  difference  formulation  we  introduce  the.  difference 
flux,  tj»(p),  as 

'l'(p)  =  $2(p)  -  ^(p),  (8) 

where  ^(p)  is  as  defined  in  Eq.  6  and  <J>^(p)  is  as  defined  in  Eq.  2. 
‘The  Boltzmann  equation  for  the  difference  flux  (using  Eqs.  2,  6,  and  8) 
is  found  to  be 

H^p)  <Kp)  =  f(p)»  ’  (9) 

where 

f(p)  =  [H1(p)  -  Hg(p)]  ^g(p).  (10) 

Assume  <J>_(p)  throughout  all  phase  space  is  known,  then  f(p)  can  be 
considered  as  a  source  which  is  non  zero  only  in  the  region  occupied 
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by  the  vehicle.  Introducing  the  adjoint  difference  flux  ty*(p),  the 
effect  of  interest  can  be  obtained  from 

X  »  /  <j>2(p)R(p)dp  -  /tf>*(p)f(p)dp,  (ll) 

where  no  approximations  have  been  introduced.  In  Eq.  11,  the  adjoint 
flux  is  the  solution  to* 

H*(p)  ij>*(p)  =  R(p)  (12) 

with  the  physical  boundary  condition  that  *J/H(p)  for  (n*fi)  >  0  goes  to 
zero  at  infinity  (or  the  physical  boundary  of  the  system). 

As  in  the  Surface  Integral  Approximation  technique,  the  Adjoint 
Difference  technique  requires  the  solution  to  the  free  field  radiation 
problem  (solution  to  Eq.  6)  and.  the  adjoint  function  which  is  the 
solution  to  Eq.  12.  Then  the  effect  of  interest  is  obtained  by  solving 
Eq.  11.  The  relative  source-target  locations  are  accounted  for  in  the 
solution  of  Eq.  11, 


*Eq.  12  differs  from  Eq.  4  only  in  the  imposed  boundary  conditions 
on  the  dependent  variable.  If  the  same  boundary  conditions  are 
imposed,  it  is  shown  in  R*>4  that  the  effect  of  interest  irom  Eqs,  7 
and  11  are  mathematically  equivalent. 


<?£ 


III.  RESULTS  FOR  NONFISSILE  PROBLEMS 

A.  Ifatroduction 

A  total  of  seven  test  problems  which  do  not  contain  fissile  materials 
have  been  considered.  The  first  three  of  these  seven  problems  can  be 
treated  as  one-dimensional;  therefore,  the  results  obtained  for  tnese 
one  dimensional  problems  from  the  Surface  Integral  and  Adjoint  Difference 
techniques  are  compared  with  more  conventional  techniques ,  The  next 
three  problems  are  geometrically  more  complex;  hence  only  the  results 
from  the  Surface  Integral  and  Adjoint  Difference  techniques  are 
presented.  The  first  six  problems  all  involve  a  target  located  in  an 
infinite  air  medium,  but  the  seventh  problem  involves  a  target  (per¬ 
turbing  region)  located  in  air  at  the  air-ground  interface.  Accord¬ 
ingly,  the  forward  flux  (the  unperturbed  flux)  for  the  first  six 
problems  was  obtained  using  the  one-dimensional  discrete  ordinates  code 

ANISN  ,  and  the  forward  flux  for  the ‘seventh  problem  was  obtained  using 

2 

the  two-dimensional  code  DOT. 

Results  obtained  from  several  techniques  will  be  presented;  therefore, 
the  next  section  io  included  for  the  purpose  of  clarifying  the  nomencla¬ 
ture  used  in  the  presentation  of  the  results.  Following  the  es uablishment 
of  the  nomenclature,  the  results  obtained  for  the  seven  problems  will 
be  presented  in  the  order  in  which  they  were  considered. 

B.  Nomenclature 

In  every  case,  the  objective  was  to  obtain  an  effect  of  interest, 

A,  as  given  by  Eq.  1,  In  some  problems,  Eq.  2  was  solved  directly 
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using  the  discrete  ordinates  code  ANISK1  and  then  solving  Eq.  1.  Results 
obtained  in  this  manner  will  be  classified  as  the  "AHISN  method." 

A  second  direct  approach  is  to  solve  Eq,  4  (the  adjoint  flux)  for 
the  overall  problem.  Then 

X  =  S(p)dp«  (13) 

P 

The  adjoint  equations  are  tolvcd  using  Monte  Carle  techniques;  therefore, 
results  obtaired  this  way  are  classified  as  the  "Adjoint  Monte  Carlo 
method."  The  ANISN  and  Adjoint  Monue  Carlo  methods  are  both  well 
established  techniques;  hence  results  obtained  this  way  are  accepted 
as  the  correct  results. 

The  use  of  Eq.  3  for  the  estimation  of  X  is  classified  as  the 
"Surface  Integral  method."  The  more  practical  but  approximate  estimate 
of  X  from  Eq,  7  is  the  "Surface  Integral  Approximate  method." 

The  use  of  Eq.  11  for  the  evaluation  of  X  is  the  "Adjoint  Difference 
method."  Finally,  the  use  of  the  <j>K(p)  of  Eq.  4  in  Eq.  11  in  place  of 
'P(p)  leads  to  the  "Adjoint.  Difference  Approximate  method." 

C.  Problems  1,  2,  and  3 

The  first  three  problems  were  selected  to  demonstrate  the  validity 
and  utility  of  the  surface  integral  and  adjoint  difference  formulations. 
Accordingly,  problems  are  considered  which  could  be  analysized  with 
the  more  standard  procedures. 

Description  of  Problem  1,  2,  and  3.  For  each  problem,  the  unper¬ 
turbed  problem  consists  of  a  point  fission-spectrum  neutron  soutcp 
located  in  an  infinite  air  medium.  Since  the  average  mean  free  path 
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of  a  fission  neutron  in  air  is  86.5  meters ,  the  unperturbed  flux  is 

assumed  spatially  constant  over  all  perturbing  regions.  The  multigroup 

transport  equations  are  used  as  he  analytic  model.  The  group 

constants  are  thore  normally  used  for  ai"  transport  problems  using 
1  2 

the  codes  ANISN  or  DOT  by  the  Neutron  Fnysics  Division  of  tie  Oak 
Ridge  National  Laboratory.  For  the  results  presented  in  this  repoit,' 
all  ANISN  calculatic.  s  were  p — -"'ormed  ucing  either  S^.  or  angular 
quadratures  for  the  flux. 

For  Problem  1,  the  perturbation  consisted  of  two  concentric 
spherical  shells  or  iron  x2  cm.  thick)  with  a  6  cm.  air  gap  between 
the  iron  shells >  The  quantity  of  interest  was  assumed  to  be  the 
response  of  a  detector  (specified  below)  located  at  the  midpoint  of 
the  air  gap  enclosed  by  the  spherical  shells.  The  detector  is  1000  m 
from  the  point  source,  fhe  geometric  details  are  presented  in  Fig.  1. 

Problem  2  was  identical  to  Problem  1  with  the  exception  that  the 
air  gap  enclosed  by  the  spherical  shells  was  filled  with  water. 

Problem  3  consists  of  a  2  cm.  thick  iron  spherical  shell  with  an 
outside  diameter  of  10  cm.  The  shell  is  filled  with  water.  A  point 
lust  flux  detector  is  located  at  the  center  of  the  spherical  shell. 

The  problem  is  to  compute  the  detector  response  at  specified  distances, 
R,  from  the  point  fission  source.  The  geometric  configuration  is 
shown  in  Fig,  2. 

Results  for  Problems  1,  2,  and  3.  The  results  fou  Problem  1  are 
presented  in  Table  1.  The  computed  detector  response  is  the  total 
fast  flux  for  the  first  13  energy  groups  (E>0,111  Mev). 


TABLE  I 


RESULTS  FOR  PROBLEM  I 


METHOD 

HISTORIES 

FLUl’RCr 

COMMENTS 

ANISN 

T.isxio'*13 

SU  A 

Adj.  Diff. 

4000 

6.95X10~13 

ANISN 
H  calc. 

Surf.  Int. 

Surf.  Int. 
Approx. 

4000 

4000 

6.98X10"13 

7.15X10"13 

The  results  for  problem  2  are  presented  in  Table  2 . 

The  results ,  total  fast  flux  for  the  first  13  groups ,  for  Problem 
3  are  presented  in  Table  3.  The  parameter  R  is  the  distance  from  the 
source  to  the  detector,  see  Fig.  2.  The  surface  over  which  the  inte¬ 
grations  were  carried  out  for  the  Surface  Integral  Approximation  method 
was  the  Fe-air  interface  surface. 

Discussion  of  Results  for  Problems  1.  2,  and  3.  From  the  results 
from  the  ANISN,  Adjoint  Difference,  and  Surface  Integral  methods  for 
problem  1  and  2  (Tables  1  and  2),  it  is  concluded  that  the  formula¬ 
tions  being  used  are  correct.  The  results  obtained  from  each  of  these 
three  methods  agree  within  statistics  as  predicted.  Foi  problem  2, 
there  is  a  slight  difference  between  results  obtained  from  the  Surface 
Integral  Approximate  method  and  the  other  methods,  but  this  difference 
was  expected  due  to  the  approximation  itself. 

Problem  3,  as  described  in  Fig.  2,  is  a  three-dimensional  problem, 
but  it  can  be  replaced  by  an  equivalent  one-dimensional  problem  by 
replacing  the  point  source  by  an  equivalent  spherical  shell  -source  of 
radius  R  whose  origin  is  at  the  center  of  the  perturbing  spherical  shell. 
This  equivalent  one-dimensional  problem  is  what  was  analysized  to  obtain 
the  ANISN  results  of  Table  3.  All  other  results  was  a  direct  analysis 
of  the  problem  as  described  in  Fig.  2,  i.e.,  the  unperturbed  flux  was 
for  the  point  source  in  an  infinite  air  medium,  and  the  adjoint  flux 
was  for  the  water  -  iron  perturbing  region  as  shown.  The  results  at  the 
different  separation  distances  R  for  the  Adjoint  Difference  and  the 
Surface  Integral  Approximate  methods  were  obtained  for  a  single  adjoint 
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SABLE  2 


RESULTS  FOR  PROBLEM  2 


Method 

Histories 

Fluence 

Comments 

MISN 

1.33  x  10~13 

Adj.  Diff, 

4,000 

1.27  x  10~13 

S4 

Surf,  Int, 

4,000 

1.20  x  10“13 

MISS  Calc. 

Surf.  Int, 
Approx, 

4,000 

1.51  x  lO-13 

■ 

ANISN 

1.64  x  10**13 

Adj .  Diff. 

50,000 

1.56  x  10“13( 0.062) 

Sl6  MISN 

Surf.  Int. 

40,000 

1.53  x  10“13 

Calc . 

Surf,  Int, 
Approx. 

40,000 

1.92  x  10~13(0.033) 

*The  numbers  in  the  parenthesis  are  the  fractional  standard 
deviation  of  the  reported  reBultB  based  on  the  Monte  Carlo  calculations. 
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Monte  Carlo  calculation  for  each  method.  The  results  for  the  Adjoint 

Monte  Carlo  method  required  a  Monte  Carlo  calculation  for  each  separation 

distance.  From  this  problem  it  was  concluded  that  the  formulations  used 

for  breaking  the  overall  problem  up  into  two  independent  separate 

problems  and  followed  by  the  combination  of  results  to  obtain  the  effect,, 
of  interest  are  correct. 

* 

D.  Problems  4,  5,  and  6 

The  unperturbed  problem  consists  of  a  unit  point  fission  source  in 
an  infinite  air  medium.  In  each  of  these  three  (unlike  the  first  three 
problems  considered)  problems,  the  effect  of  interest  is  dependent  on 
the  orientation  of  the  perturbing  medium  (vehicle)  relative  to  the 
source.  Problems  4  and  5  consists  of  the  same  vehicle  with  different 
detectors.  In  the  analysis  which  were  carried  out  with  Problems  4  and 
5,  it  was  found  that  a  relatively  large  ©mount  of  computer  time  was 
required  due  to  fche  use  of  the  generalized  geometry  package  in  MORSE; 
therefore,  Problem  6  was  introduced  which  has  the  same  general  features 
of  Problems  4  and  5  and  does  not  require  the  more  generalized  geometry 
routines.  Therefore,  statistically  more  precise  results  were  obtained 
•Po**  Pj*q!xLs21  &  • 

Description  of  Problems  4,  5,  and  6,  In  Problem  4,  the  perturbing 
medium  was  taken  to  be  the  vehicle,  made  of  iron,  shown  in  Fig.  3.  The 
effect  of  interest  was  the  fast  flux  for  the  point  detector  located  as 
shown  in  Fig.  3.  Problem  5  differs  only  from  Problem  4  through  the  fast 
flux  detector.  In  particular.  Problem  5  consists  of  a  six  inch  diameter 
spherical  detector  centered  at  the  same  point  as  the  point  detector  of 
Problem  4. 
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Three  different  orientations  of  the  vehicle  (called  rocket  from 
henceforth)  are  considered  for  Problems  4  and  5.  The  *e  are: 

1.  The  rocket  is  moving  around  the  source  at  a  distance  R 

away  such  that  the  z»axis  is  perpendicular  to  a  straight 
line  connecting  the  source  point  and  the  detection  point. 

2.  The  nose  of  the  rocket  is  moving  away  from  the  source. 

3.  The  nose  of  the  rocket  is  moving  toward  the  source. 

Problem  6  consists  ex'  the  iron  parallelpiped  shown  in  Fig.  It.4* 

The  effect  of  interest  is  the  Bame  as  in  the  previous  problems  (the  fast 
flux) .  The  fast  flux  is  desired  for  three  orientations  of  the  parallel- 
piped  at  each  specified  distance  between  the  source  and  the  detector. 

These  orientations  are  illustrated  by  the  Roman  numerals  in  Fig.  4. 

The  numerals  represent  the  source  direction  relative  to  the  parallelpiped. 

Results  for  Problems  it,  5,  and  6.  The  computed  fast  flux  for 
Problems  4  and  5  and  presented  in  Tables  4  and  5  for  each  of  the 
indicated  rocket  orientations.  The  fractional  standard  deviations  in 
all  of  the  results  presented  in  Tables  4  and  5  are  on  the  order  of  25#. 

The  fast  flux  multiplied  by  4ttR2  (R  is  the  distance  in  meters  from 
the  source  to  the  detector)  obtained  for  Problem  6  are  presented  in 
Table  6.  The  fractional  standard  deviation  is  less  than  10#  for  all 
casej  presented  for  Problem  6  as  indicated  in  Table  6*. 

^Problems  4  and  5  were  found  to  consume  a  large  amount  of  computer 
time  due  to  the  generalized  geometry  package;  therefore,  we  choose  the 
parallelpiped  to  continue  the  source-vehicle  orientation  study.  This 
geometry  requires  relatively  little  computer  time. 
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TABLE  4 

RESULTS  FOR  PROBLEM  4a 


Distance  from 

Source  Me  tors _  Ad,1.  Dtff, _ Surf,  Int.  Approx. 


Around  Scurce 
2C0 
1*00 
600 
800 
1000 


3.61*1  x 
5.128  x 
1.026  x 

2.25I*  X 
5.035  x 


10 

10' 

10' 

10' 

10 


■10 

-11 

•11 

-12 

-13 


3.1*71*  x  10"?-° 
5.565  x  10““ 
1.100  x  10,1 
2.39.0  x  10",, 
5.560  x  10'  J 


Away  From  Source 
200m 
1*00 
600 
800 
1000 


2.937  x 
4.1*30  x 
9.067  x 
2. Oil*  x 
4.527  x 


10 

10' 

10' 

10' 

10' 


-10 

-11 

-12 

-12 

-13 


3.137  x  10'^ 
5.192  x  10,, 
1.046  x  10'  ,r 
2.315  x  10 

5.421  x  10  J 


Toward  Source 
200m 

400 

600 

800 

1000 


3.251  x 
4.8l6  x 
9.829  x 
2.182  x 
4.904  x 


10 

10 

10 

10 

10 


-10 

-11 

-12 

-12 

-13 


4.389  x 
6.508  x 
1.296  x 
2.829  x 
6.5l6  x 


10 

10 

10 

10 

10 


-10 

-11 

-11 

-12 

'-13 


£L 

ANISN  calculations  were  S^. 
Based  on  4,000  histories. 


Estimates  for  Adjoint  flux  were 
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TABLE  5 

RESULTS  FOR  PROBLEM  5& 


Distance  from 

Souice  Meters 

AdJ .  Diff. 

Surf.  Int.  Approx. 

i 


Around  Source 
200 
4oo 
600 
800 
1000 


4.390  x 
6,84c  x 
1.433  x 
3.267  x 
7.830  x 


10 

10' 

10' 

10 

10' 


.-10 


£ 

-12 

-13 


3.448  x  10~|? 
5.568  x  10“:r 
1.172  x  lO^ 
2.609  X  io”r* 

6.459  X  10  J 


Away  From  Source 


200m 

3.513 

400 

6.124 

600 

1.329 

800 

3.078 

1000 

7.435 

x  10"|°  2.920  x  10“|° 

x  io~rr  4.750  x  io"rr 

X  10“rr  1.031  X  10 

x  io"r^  2.4o4  x  io“r^ 

x  10"15  5.830  x  10'  * 


Toward  Source 

200m 

3.858 

400 

6.580 

600 

1.420 

800 

3.281 

1000 

7.916 

x  10"|°  3.976  x  10“|° 
X  io”rr  6.463  X  io“r: 
x  io“rf  1.351  x  io“|i 
x  io~:r  3.071  x  io“r^ 
x  10  J  7.327  x  10  5 


aAll  results  are  based  on  5»000  adjuncton  histories  and  S  ^  angular^ 
fluxes.  The  tabulated  values  are  spatially  averaged  fast  flux  fluences. 
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TABLE  6 

RESULTS  FOR  PROBLEM  6 
(4ttR2  «  Fast  Flux) 


R  (Meters) 


Ad,) .  Mff.  Surf.  Int.  Approx. 


Orientation  X 
200 
4oo 
600  , 

800 

1000 

Orientation  II 
200 
4oo 
600 
800 
1000 

Orientation  III 
200 
1*00 
600 
800 
1000 


2.268  (0.05l)a 
1.336  (0.034) 
0.619  (0.029) 
0.250  (0.027) 
0.0931(0.026) 


1.772  (0.075) 
1.156  (o.o46) 
0.558  (0.038) 
0.229  (0.035) 
0.0864(0.033) 


1.563  (0.088) 
1,078  (0.053) 
0.525  (0.045) 
0.216  (o.G4l) 
0.0815(0.039) 


2.268  (0.056) 
1.369  (0.058) 
O.636  (0.059 
0.255  (0.059) 
0.0951(0.060) 


1.819  (0.058) 
1.171  (0.058) 
0.  61  (0.058) 
0  228  (0.059) 
0.0857(0.059) 


1.590  (0.059) 
1.077  (0.060) 
0.524  (0.061) 
0.214  (0.062) 
0.0810(0.063) 


aFractional  standard  deviations  in  parentheses. 


The  calculational  procedure  followed  in  Problems  4  and  5  differed 
from  that  pursued  in  Problem  6,  In  particular,  the  procedure  used  in 
Problems  4  and  5  was  the  three  step  procedure  outlined  in  Sec.  Hof  this 
report.  The  technique  employed  in  Problem  6  was  to  evaluate  Eq.  7  (for 
the  Surface  Integral  Approximation)  or  Eq.  11  (for  the  Adjoint  Difference) 
directly  in  the  Monte  Carlo  code  MOPSE. 

Discussion  of  Results  for  Problems  4,  5,  and  6.  From  the  results 
obtained  for  Problems  4  and  5,  it  is  concluded  that  the  three-step 
procedure  outlined  earlier  in  this  report  is  valid.  In  particular,  we 
can  obtain  the  effect  of  interest  within  the  target  for  different  source- 
target  orientations  by  using  a  single  forward  flux  calculation  for  the 
source  in  air  problem  and  a  single  adjoint  flux  calculation  for  the 
target  in  air  problem.  However,  it  was  also  found  for  a  generalized 
application  of  this  three-step  technique,  further  development  iB 
required  to  properly  interface  the  results  obtained  from  the  two  indepen¬ 
dent  calculations.  The  problems  are  introduced  by  (a)  the  use  of 
different  reference  coordinate  systems  in  the  two  calculations  and 
(b)  coupling  the  results  obtained  from  deterministic  calculational 

fonhnn  rmoa 
— u — -  . 

The  results  for  Problem  6  were  obtained  by  (a)  evaluate  the  "response 
function"  using  the  forward  flux  and  the  coordinate  system  used  in  the 
forward  flux  calculation,  (b)  transform  the  response  function  to  the 
coordinate  system  to  be  used  in  the  adjoint  (Monte  Carlo)  flux  calcula¬ 
tions,  and  (c)  evaluate  the  effect  of  interest  directly  in  the  Monte  Carlo 
code  using  the  transformed  response  function.  The  disadvantage  with 
this  procedure  is  that  an  adjoint  calculation  is  required  for  each 


target-vehicle  orientation. #  This  approach  was  adopted  to  (a)  circumvent 
the  problems  encountered  in  the  Problems  4  and  5»  and  (b )  to  establish 
a  reference  calculations!  technique  for  the  more  generalized  techniques 
which  are  to  be  developed  to  utilize  the  three  step  calculational 
techniques. 


E.  Problem  7 

Description.  The  unperturbed  problem  is  thau  of  calculating  the 
angular  fluxes  due  to  a  unit  point  fission  source  located  50  feet  above 
an  air-ground  interface.  This  is  a  two-dimensional  (r-z  cylindrical 
geometry)  problem.  This  problem  was  solved  by  E.  A.  Straker’*  with  the 
discrete  ordinates  code  DOT.  The  unperturbed  fluxes  used  in  this 
problem  were  taken  from  Straker's  work. 

The  perturbation  to  this  system  is  a  two-foot  right  concrete 
cylinder.  At  the  center  of  this  cylinder  is  a  three  inch  cubic  detector. 
The  detector  measures  the  average  fast  fluence  (E>.1  Mev)  in  the  detector 
volume.  The  geometric  configuration  is  illustrated  in  Figure  5.  The 
concrete  (with  rebar)  and  ground  (8.5  weight  percent  water)  compositions 
are  shown  in  Table  7. 

Scoring  Surface  for  the  Conventional  Method.  In  previous  problems, 
the  surface  used  in  the  conventional  method  was  the  surface  of  the 


#This  technique  is  still  preferable  over  the  more  direct  approach  of 
explicitly  accounting  for  the  source  and  target  in  a  single  calculation. 

The  difficulty  with  the  single  calculation  is  the  normally  large  separation 
between  the  target  and  source  which  is  a  deep  penetration  geometrically 
complex  problem. 


I 
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conventional  method 


TABLE  T 


MATERIAL  COMPOSITIONS 


(Number  Density  x  10  ) 

Element 

Concrete  &  Rebar 
(2.62  gm/cm^) 

Ground 
(1.6  gm/cm^) 

0 

4,084 

E-2 

3.036 

E-2 

Si 

1.332 

E-2 

9.990 

E-3 

A1 

2.738 

E~3 

2.04Q 

E-3 

Fe 

5.269 

E-3 

2.400 

£-4 

Mg 

1.620  E-4 

1.200 

E-4 

Ca 

2.426 

E-3 

1.800 

E-3 

Na 

1.071 

E-3 

8.000 

E-4 

K 

8.280 

e-4 

6.200 

E-4 

H  ■ 

1.065 

E-2 

7.600 

E-3 

Ti 

2.600 

E-5 

2.000 

E-5 

Mn 

•  7.500 

E-5 

2.000 

E-5 

S 

8.400 

E-5 

6.000 

E-5 

C 

1.310 

E-4 

-  - 

P 

5.000 

E-6 

•urn  mm  ■ 

-  - 

Hi 

1.300 

e-4 

-  -  . 

— 

Cr 

1.200 

E-5 

-  - 

Mo 

1.800 

E-5 

- - 

■■  - 

Cu 

1.500 

E-5 

—  ~  - 

-  - 

■26- 


I"  I 

t  •( 


k> 


% . 
f" ! 

V. 

I  * 
\\ 


*  t 

NT  J 

H 


si 

I 

b 

iv 

& 

fe 

fl 

li 

U 


i 

i « 


r, 

li 


perturbing  material.  In  this  problem,  however,  the  inward  directed 
perturbed  flux  will  be  considerably  less  than  the  inward  directed  unper¬ 
turbed  flux  at  the  surface  of  the  cylinder  in  contact  with  the  air-ground 
interface.  Therefore,  the  surface  sketched  in  Figure  4,  which  extends 
in  the  vertical  direction  to  minus  infinity,  was  used. 

The  forward  and  adjoint  fluxes  decrease  rapidly  with  distance  below 
the  air-ground  interface,  (The  mean  free  path  in  ground  is  about  two 
centimeters.)  Therefore,  the  contribution  to  the  effect  of  interest  in 
the  conventional  method  of  the  surface  below  the  interface  was  assumed 
insignificant. 

Results  for  Problem  2.  The  unperturbed,  DOT,  calculations  were  Sg 
(40  discrete  directions).  The  cross  sections  used  in  DOT  and  MORSE  were 
22  group,  Pg.  The  adjoint  Monte  Carlo  calculations  involved  following 
20,000  histories.  Since  the  perturbation  greatly  reduced  the  fast  fluence, 
the  change  in  the  effect  of  interest  is  reported  in  Table  8,  i.e.,  for 
the  adjoint  difference, 

AX  «  /  **(p)  f(p)dp.  '  ’  (14) 


and  for  the  conventional, 

ua  »  /  R(y)  <j*p(p)dp  +  /  (n  8  n).  $p(p)  rw(p)dp.  (15) 

Det  Surface  in  * 

Discussion  of  Results  for  Problem  7.  In  Table  3,  the  obvious 
difference  between  the  surface  integral  approximation  and  the  adjoint 
difference  method  is  the  statistical  uncertainty.  In  particular,  the 
fractional  standard  deviations  for  the  surface  integral  approximation 
are  observed  to  be  very  low.  This  effect  can  be  accounted  for  by  noting 
that  Eg,.  (15)  consists  of  adding  a  term  calculated  by  Monte  Carlo 
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techniques  to  a  term  calculated  "by  discrete  ordinates  techniques.  For 
large  perturbations »  as  was  the  case  for  this  problem,  the  term  calculated 
by  Monte  Carlo  is  small,  Therefore,  even  a  relatively  large  fractional 
standard  deviation  vas  in  the  Monte  Carlo  calculation  will  appear  small 
overall.  This  effect  is  not  observed  in  the  adjoint  difference  method 
because  the  statistical  uncertainty  of  the  quantity  calculated  with  Eq. 
(lH)  is  a  measure  of  the  reproducibility  of  the  Monte  Carlo  calculation 
alone , 

The  results  from  the  two  calculational  techniques  (surface  inte¬ 
gral  and  adjoint  difference)  agree  with  each  other  well.  This  was 
expected  since  little  doubt  exists  as  to  how  to  choose  the  surface  for 
the  surface  integral  technique,  hut  for  some  problems  of  this  type, 

(a  vehicle  located  at  an  air-ground  interface  for  example)  it  may  be 
difficult  to  choose  the  surface  apriori.  In  such  a  situation,  the 
adjoint  difference  technique  would  be  preferable  over  the  surface 
integral  technique, 
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TABLE  8 


RESULTS  FOR  PROBLEM  ?& 


R 

Adjoint  Difference 

Conventional 

(Meters ) 

Method 

Method 

300 

0.893  (0.078)b 

0.882  (0.008) 

600 

0.21+9  (0.065) 

0.253  (0.008) 

900 

0.0603  (0.061) 

0.0620  (0.008) 

1200 

0.0135  (0,057) 

0.0139  (0.008) 

1500 

0.00315(0,056) 

0.00325(0.008) 

aAll  results  are  (UrrR^  *  change  in  fast  fluence )  . 
Fractional  standard  deviation  in  parentheses. 
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IV.  RESULTS  FOR  FISSILE  PROBLEMS 


The  re a \ilts  for  a  single  fissile  problem  are  included  in  this  report. 
The  objective  is  to  demonstrate  the  applicability  of  the 'adjoint  mode  of 
the  MORSE  code  to  problems  of  the  type  of  interest  which  contain  fissile 
material.  The  problem  chosen  is  the  "Calculation  No.  2"  problem 
presented  in  Ref.  6.  The  results  will  be  compared  with  ANISN,  and  the 
ANTE-2  and  DTF  results  presented  in  Ref.  6. 

The  problem  is  to  compute  the  number  of  fissions  per  kilogram  of 
U^J^  in  the  target  shown  in  Fig.  6.  The  spherical  target  was  subjected 
to  a  parallel  incident  neutron  beam  as  shown.  The  incident  neutron  beam 
energy  range  from  lb  Mev  to  thermal  was  broken  up  into  a  25  group  structure 
as  given  in  Table  9.  The  fission  cross  section  given  in  Table  9  was  used 
as  the  response  function.  The  results  are  presented  in  Fig.  7  as  a  func¬ 
tion  of  the  incident,  neutron  energy  over  the  energy  range  covered  from 
group  1  through  group  18. 

The  discrepancies  between  the  ANTE-2  and  DTF  results  were  attributed 
to  the  different  cross  sections  used.’’  The  same  cross  section  set  was 
used  in  MORSE  and  ANISN  (but  different  than  those  used  in  Ref.  6). 

Therefore  j  we  conclude  that  the  agreeiueut  between  all  focus  calculations.! 
tools  are  as  good  aa  can  be  expected.  From  these  results,  it  is 
concluded  that  the  adjoint  mode  of  MORSE  can  be  applied  to  problems 
containing  fissile  materials.  The  cnlculational  procedure  employed  in 
this  sample  problem  was  basically  the  surface  integral  approximation; 
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TABLE  9 

ENERGY  BAND  STRUCTURE  FOR  FISSILE  SAMPLE  PROBLEM 


Group 

Energy 

1 

12  -  l4.  Mev 

1.92967 

2 

8.3  -  12 

1.79722 

3 

5.3  -  8.3 

1.36285 

4 

3.4  -  5.3 

1.25280 

5 

2.2  -  3.4 

1.30383 

6 

1.4  -  2.2 

1.24789 

7 

0.9  -  1.4 

1.18587 

8 

.58  -  0.9 

1.15389 

9 

370  -  580  kev 

1.17654 

10 

240  -  370 

1.26027 

11 

150  -  240 

1.39574 

12 

100  -  150 

1.57489 

13 

31  -  100 

2.25959 

l4 

10  -  31 

3.15976 

15 

3.16  -  10 

4.20512 

16 

1.0  -  3.16 

6.37329 

17 

.316  -  1.0 

11 . 4i85 

18 

100  -  316  ev 

'20.3925 

19 

31.6  -  100 

33.9860 

20 

10  -  31.6 

57.2095 

21 

3.16  -  10 

38.5393 

22 

1.0  -  3.16 

28.3100 

23 

.316  -  1.0 

84.6518 

24 

.076  -  .316 

241.256 

25 

.015  -  .026 

577.347 

32. 


(upper set)  INCIDENT  ENERGY  »v> 


therefore,  we  have  established  that  the  techniques  applied  to  the  non- 
fissile  problems  can  be  extended  to  fissile  problems.  The  specific 
problems  encountered  will  be  reported  at  a  later  date  along  with  specific 


V.  SUMMARY 


In  this  report}  a  brief  theoretical  development  was  presented  for 
the  adjoint  difference  and  surface  integral  techniques  for  solving  a 
certain  class  of  deep  penetration  geometrically  complex  radiation 
transport  problem.  The  reader  was  referred  to  Ref.  4  for  a  more 
thorough  theoretical  development.  In  Ref.  4,  a  selected  amount  of 
calculated  results  for  targets  were  presented  to  demonstrate  the 
accuracy  of  the  theoretical  developments.  In  this  report,  all  results 
obtained  in  this  study  for  targets  which  do  not  contain  fissile 
materials  are  presented.  Furthermore,  some  preliminary  results  are 
presented  for  a  target  which  does  contain  fissile  materials  are 
presented  and  compared  with  results  reported  in  the  literature. 

From  the  overall  results  presented  in  this  report,  it  is  obvious 
that  the  adjcint  difference  and  surface  integral  methods  both  lead  to 
acceptable  results.  In  general,  it  appears  (based  on  results  from 

I 

the  non- fissile  problems)  that  the  adjoint  difference  method  is  advanta¬ 
geous  for  those  problems  in  which  the  target  introduces  a  small  pertur¬ 
bation  and  the  surface  integral  approximation  is  advantageous  (provided 
an  adequav 2  "scoring"  surface  can  be  identified  apriori)  for  those 
problems  in  which  the  target  introduces  a  large  perturbation.  These 
conclusions  cannot  be  carried  over  to  problems  containing  fissile 
materials  at  this  time  (they  will  be  explored  via  examples  and  reported 
at  a  later  date). 

At  this  time,  the  utility  of  the  three  steps  presented  in  the  intro¬ 
duction  has  not  been  fully  realized  due  to  the  complexity  of  Step  3  for 
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a  target.  This  phase  of  the  overall  problem  will  he  explored  further 
and  reported  at  a  later  date. 
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